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I.  INTRODUCTION 

In  earlier  reports  (1,  2),  a  method  for  calculating  the  thermal 
resoonse  due  to  intensive  heating  was  described.  The  present  report 
describes  extensions  and  applications  of  the  method  to  other  laser 
heating  problems.  In  the  second  chapter,  a  general  criterion  for  ore 
dimensional  heating  under  the  assumption  that  the  melt  is  instantaneously 
removed  will  be  developed.  The  predicticnsof  a  simple  analytical 
solution  for  radial  conduction  are  combined  with  these  results  in  the 
cmrd  chapter  to  produce  a  method  for  estimating  from  a  single  curve 
the  melt  through  time  for  any  arbitrary  combination  of  flux,  spot  size, 
thickness  and  material.  In  the  fourtr.  chapter, melt  through  time*  is  computed 
for  cne  case  of  complete  retention  of  the  melt  until  vaporization  occurs> 
borne  predictions  for  neating  due  to  a  large  single  pulsje  - 

will  be  presented  in  the  fiftn  chapter  and  compared  with  the  results  of 
experiments  featuring  large  single  pulses. 
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II.  A  CRITERION  FOR  ONE  DIMENSIONAL  MELTING -OF  SOLIDS 

The  question  of  the  flu<  (or  power  level)  required  to  generate  one- 
dimensional  (axial)  melting  is  particularly  significant  in  the  design 
and  interpretation  of  experiments  on  laser  effects.  Radial  heat  conduc¬ 
tion  has  been  found  to  make  the  interpretation  of  some  experimental  re¬ 
sults  difficult.  The  foiling  analysis  leads  to  a  criterion  for  one- 
dimensionality  which,  it  is  hoped,  will  prove  useful  in  the  design  of 
experiments. 

Ue  assume  that  a  melting  front  moves  at  a  constant  speed  V  through 
a  slab  of  material  having  thermal  properties  K,  p  and  Cp,  and  a  melting 
temperature  Tm.  Material  at  distances  greater  than  some  length  i  will 
essentially  be  at  an  ambient  temperature  TQ,  i.e., 

•<5>-V<VV  *  *  0) 

We  assume  that  the  heat  flow  can  be  regarded  as  one  dimensional  (axial) 

if  the  depth  of  heating  is  small  compared  to  the  width,  or  if 

6/o<N  (2) 

where  a  is  seme  characteristic  dimension  of  the  heated  region,  and  N  is 

assumed  to  be  independent  of  flux  intensity  and  thermal  properties.  In 

essence,  we  assume  one-dimensionality  if  temperature  changes  take  place 

over  distances  much  less  than  the  dimension  of  the  heated  area.  If  the 

heat  flow  is  one  dimensional ,  the  heat  conduction  equation  reduces  to 

32T  ,  3T 
3X2  3t 

where  the  diffusivity,  x  is  defined  by 

.  K 


(3) 


pcr 


(4) 
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s 

> 
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and  x  is  taken  to  he  normal  to  the  surfaces  of  the  slab.  In  terms  of  a 


coordinate  system  moving  with  the  front,  i.e.. 


Equation  (2)  becomes 


2  *  x-Vt 


K  2fl.  .  v  II 

3Z2  32 


which  has  the  solution 


From  1, 


T  *  T0  ♦  (Tm-T0)e~VZ/* 


6V/k  *  in  (1/c) 


(5) 

(6) 

(7) 

(8) 


One  dimensional  axial  heat  flow  may  also  exist  in  a  sheet  so  thin 
as  to  be  of  nearly  uniform  temperature,  i.e.,  if  the  characteristic  dimen¬ 
sion,  o,  of  the  heated  region  is  large  compared  with  t,  the  thickness,  or 

i/o<M  (9) 

M  is  a  dimensionless  number,  assumed  to  be  independent  of  material  prof 
erties.  More  generally  then,  we  expect  that  the  heat  flow  will  be  pre¬ 
dominantly  axial  (or  one  dimensional),  if  either  of  inequalities  1  and  9 
are  satisfied,  or  if  a  new  inequality  combining  the  two  is  satisfied,  i.e. 


if 

16/0*<M*  (10) 

where  M*  is  assumed  to  be  independent  of  flux  intensity  aid  thermal 
properties. 

The  steady  speed,  V,  at  which  the  liquid-solid  interface  moves  may 
be  readily  determined  from  a  heat  balance.  In  any  time  At,  a  thickness 
A2  of  material  at  a  temperature  arbitrarily  close  (c)  to  the  ambient  temp¬ 
erature  Tq  enters  the  heated  layer.  During  that  time  increment,  a  like 
amount  of  material  is  melted.  Thus,  for  any  unit  area,  a  thickness  az 
is  raised  from  the  ambient  to  melting  temperature,  and  undergoes  phase 
transformation. 
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The  heat  required  is  supplied  by  the  absorbed  intensity  Ip3,  and  a  simple 
heat  balance  on  a  unit  area  subject  to  the  peak  intensity  yields 

Ipaat  «  [Cp(Tm  -T0Hm]paz  (11) 

where  Lq  is  the  heat  of  fusion.  Thus  if  Jpa  (the  peak  absorbed  energy 
per  unit  area  per  unit  time)  Is  constant,  the  speed  at  which  the  liquid- 
solid  interface  penetrates  the  solid  portion  at  the  point  of  maximum 
intensity  is  given  by 


AZ 


Ipa 


“  P[Cp(VT0).lJ 

Substituting  this  result,  and  Equation  8  into  inequality  10,  we  find 


02) 


t«Cp 


»2I 


CCp(Tm-T0)4Lm]  < 


M* 


pa 


in(l/e) 


(13) 


If  we  assume  the  beam  profile  to  be  Gaussian  and  axi -symmetric,  the 
distribution  cf  absorbed  intensity  Is 


I(r)  ■  V 

the  total  power  absorbed  is  given  by 

Pa  -  2»o2  I 

and  inequality  13  may  be  written  as 

^JL  >  m* 


-r2/(2o2) 


pa 


(14) 


(15) 


06) 


This  furnisr.es  a  theoretical  basis  for  the  criterion  for  one  dimensional 
axial  heat  flow  deduced  in  earlier  work  [?3  from  the  results  of  computer 
experiments  on  magnesium  where  m*  was  found  to  be  40  KW/cm. 

Equation  13  suggests,  however,  the  existence  of  a  single  constant,  a 
dimensionless  power  p^r  unit  thickness,  applicable  to  any  material.  Let 

o  —  Pa 


ana 


**n[L^Cp(Tm-T0)] 
p  *  a  2tttn( 1/e) 


M* 


0?) 


(10) 
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Then  the  criterion  for  one  dimensional  celling  becomes 

pt.  >  <’««’  09) 

Tne  assuaption  of  some  profile  other  than  Gaussian  (e.a.t  a  “flat-top"), 

could  affect  at  most  the  numerical  value  of  the  ditoensionless  Pti*. 


Equation  13  *ay  be  written  in  another  fora,  through  the  use  of  the 
tine  tj  required  for  complete  melt  through  under  the  one-dieensional 
condition.  The  required  tine  is  obtained  froa  a  one-disensfonal  heat 
balance  as  / 


f  ,  rCptTjj-TgjH^lpt 

1  - nz - 


Substituting  into  13  leaves  a  diaensionless  parameter  known  as  a 


Fourier  mmber. 


4  <  H*/in{l/c) 


Ke  have  seen  that  the  conditions  necessary  for  one-dieensional  Belting 
are  oet  if  a  Fourier  number,  defined  through 

F  «  Kt}/a2  (22) 

is  less  than  some  critical  value,  or  if  the  diirer.sionless  power  per  unit 
thickness,  defined  through  Equation  17  and  related  to  the  Fourier  number 
through 


is  greater  than  some  critical  value.  It  is  at  this  point  unproven  that 

the  critical  value  of  these  dissension! ess  parameters  are  independent  of 

material  properties,  although  the  results  given  e»»*lier  for  magnesium 

are  a  satisfactory  indication  that  the  critical  value  of  the  parameter  P, 

w 

is  independent  of  the  flux  intensity,  dimensions  of  the  heated  area,  and 
sample  thickness. 


!n  order  to  »sc«rUifi  if  the  critic*!  value  of'P^,  the  dinessionless 
abscrt^d  power  per  unit  thickness,  is  iiitepen&t/i t  of  sas^rial  paraaaiters* 
calculations  using  the  previously  described  fl&tcticz’i  estfedd  *ere  under- 
taken.  Katerial  properties  as  given  in  Table  »  were  a$»aed,  with  re¬ 
sults  as  given  in  Table  11.  In  all  cases,  a  disk  2.5  c*  dias&ter  .1  c* 
tfcscfc  was  divided  into  20  annular  rings  an*?  10  layers.  The  parasK'^r  « 
describing  Use  gauss i an  tees*  ;*?*  take*  to  be  .26  as.  The  respite  ere 
plotted  in  dimensionless  fore  as  figure  !• ,  and  suggest  that  for  ab¬ 
sorbed  dimensionless  power  per  unit  thickness  of  70  or  greater,  the 
zxltfsg  tfae  will  he  within  ICS  of  the  value  calculated  fre*  tie  cse- 
dimensional  heat  balance  vtgualioo  %Z). 
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THERMAL  PROPERTIES  OF  SEVERAL  MATERIALS 


gm 

3.6 


2.73 

1.77 


7.9 

4.43 


m 

Joy  I e/gm 
1070 

.375 

338 


2  SO 
390 


C, 


Joule/(cm  sec  °K) 
.104 
1.84 
.96 


'P 
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.24 

.145 


gm°K 

.855 


1.0 

1.04 


’m 

#K 

2313 

911 

905 


.42 

.77 


1700 

1900 


Table  II 


COMPUTED  MELT  THROUGH  TIMES  FOR  VARIOUS  LEVELS  OF  ABSORBED  FLUX 


Material 

r  iua  . 

l  Watts/or 

tm-sec 

t,-sec 

A 1 2^3 

1000 

1.4197 

1.084 

2000 

.6226 

.5418 

3000 

.3985 

.3612 

5000 

.232 

.2167 

10000 

.1140 

.1084 

30000 

.0377 

.03612 

Aluminum 

2000 

.3131 

.1371 

3000 

.1518 

.09137 

5000 

.0733 

.05482 

10000 

.0318 

.02741 

Magnesium 

1000 

.4128 

.1712 

1500 

.1962 

.1141 

2000 

.1266 

.0856 

3000 

.0737 

.05706 

5000 

.0402 

.03424 

10000 

.0188 

.0171 

30000 

.0061 

.00571 

Stainless  Steel 

1000 

1.1068 

.6936 

3000 

.270 

.2312 

5000 

.1532 

.1387 

10000 

.0742 

.06936 

30000 

.0244 

.02312 

Titanium 

300 

9.098 

2.3  73 

'  1000 

I.GC  • 

.7185 

3000 

.2/1 

.2395 

5000 

.  ,57 

.1437 

10000 

.078 

.07185 

30000 

.0255 

.02395 

ALUMINUM 


-DIMENSIONLESS  POWER  PER  UNIT  THICKNESS  (Absorbed) 


III.  MELTING  TIMES  AT  LOW  POWER 


$ 


Up  to  this  point,  only  the  times  necessary  for  complete  melt 
through  have  been  considered  for  cases  where  the  incident  flux  is 
relatively  high,  i.e.  the  heat  transfer  primarily  axial,  further 
insight  may  be  gained  through  consideration  of  the  other  extreme, 
i.e.  assuming  that  the  heat  flow  is  primarily  radial. 

Carslaw  and  Jaeger^have  given  the  temperatures  distribution 
in  a  solid  due  to  the  application  of  a  flux  Q  per  unit  area  applied 
over  the  interior  of  a  cylindrical  hole  in  a  infinite  domain.  The 
temperature  is 

T(r,t)-T  «  -  ^  r  (l-e”*u  *) 

0  wk  0 

where  R  is  the  radius  of  the  hole 
«  is  the  diffusivity 
k  the  conductivity 

J  and  Y  are  bessel  functions  of  first  and  second  kind 
and  TQ  is  the  initial  temperature 

3- 

For  large  values  of  time,  ar.  approximation  is 


where  C  =  3  1.781.  This  approximation  has  been  used  to  predict  melting  times 

for  thin  sheets  exposed  to  small  lasers.**  Retaining  only  the  first  term  leaves 


J0(ur)Y1(uR)-Yn(ur)J1(uR)  du 
J^uR)  +  Y^(uR)  u2 


T-T 


e 


.01 

2k 


;.n 


(26) 


IP 


-  -  - - - - 


I 

A 

y 


a 

! 

I 


as  an  approximate  temperature  in  a  thin  sheet.  In  order  to  apply  this 
solution  to  the  problem  of  laser  melting,  we  assume  a  uniform  flux  of 
intensely  FQ  to  be  applied  over  a  radius  R  on  a  sheet  of  thickness  t.  For 
a  small  t,  we  assume  the  temperature  under  the  beam  (i.e.  r<R)  to  be 
uniform  and  equal  to  the  temperature  at  R.  For  moderate  values  of  flux, 
we  assume  a  quasi  steady  temperature  under  the  b^am,  viz.,  t^t  the  transient 
term  may  be  lisregarde-’  Hence,  a  heat  balance  on  the  mass  under  the 


spot  read* 


F0  *  kR  *  Q  *  2lR  *  1 


where  Q  is  the  radial  flux,  as  in  Equations  24  and  26.  Using  Equation  26  we  then 
have  as  a  temperature  distribution 


R4F„  4ct 


T(r,t)  *  Tn  ♦  — a  tn 


4kl  CR2 


on  C<r<R 


(28a) 


C  d2  d**- 

T(r,t)  *  T  +  —  in  — -  on  r>R 

0  4kt  Cr2 


(28b) 


It  should  be  noted  that  certain  assumptions  made  in  the  development 
of  this  result  are  not  completely  consistent.  The  assumption  of  long 
time  and  low  intensity  are  not  compatible  with  the  neglect  of  losses 
due  to  convection  and  radiation,  hence  the  temperatures  predicted  bv  28b  are 
at  best  upper  bounds,  and  melting  times  which  might  be  predicted  are 
lower  bounds. 

Such  melting  times  may  be  estimated  from  Eouation  28a  as 


t  *  £SL»  GXD  (T  -T  ) 

P  FR2  0 
0 


Thi;  melting  time  can  be  compared  with  that  predicted  for  one  dimen¬ 
sional  axial  flow  with  instantaneous  removal. 


hence 


^  +c  (t  -tjI 

1  F  (  m  p'  m  o  5 
ro 

S»=P  5- exp  f _ ^ _ l 

‘i  * 4*  L  r. .  ^  i ) 


p  11  + 


'  L  CP!VT0>J 


where 


P  «  »*  fo  _  ( 

*  »xp[L  K  (T  -T  )] 
n  p  m  o 

The  dimensionless  power  per  unit  thickness  is  again  seen  to  characterize 
the  melting  time  ratio.  For  small  P^,  6  becomes  large.  For  large  Pt, 


but  the  assumption  of  one  dimensional  radial  flow  is  inapplicable  in  this  reolme.  J 
Equation  31  has  a  mimimum  at  I 


C  (T  -T  ) 
p  m  o 


C  (T  -T  ) 
p  m  o 


and  should  not  bo  used  for  values  of  P  above  that  given  by  Equation  34, 
The  previous  results  were  deduced  under  the  assumption  of  a 
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"flat-top"  beam  profile.  For  a  Gaussian  beam,  the  assumption  that 
the  material  under  the  beam  (R  =  2a)  melts  uniformly  is  somewhat 
questionable,  but  will  enable  comparison  with  the  results  obtained  by 
the  numerical  method.  Equation  27  then  reads 


1  F  *R2  =  Q2*Rt 

2  o 


and  Equation  29  becomes 


t  *  £a_  exp  f  (T -TJ 
Bi  <  /  f  02  m  0 


normalizing  by  the  one  dimensional  melt  through  time  (Equation  20)  and 
introducing  the  dimensionless  power  per  unit  thickness  for  a  gaussian 


p  »  2w~*pa 

*  XnVC  (T-TH 
m  p  m  o 


leaves 


_ i _ ] 

‘i  2*  Kr, ,  L-  i  \ 
.  yvv  . 


A  minimum  again  occurs  at  the  value  of  P£  which  make  the  exponent  of  Equation 
39  unity,  and  leaves 


°min 


atPt  = 


C  (T  -T  ) 
p  m  o' 


C  (T  -T  ) 
p  m  o 


The  parameter  1  +  L^/Cpd^-^)  ranges  from  1.613  in  the  case  of  aluminum 

to  1,31b  in  the  case  of  Titanium,  v/i tii  intermediate  values  of  1.537  for 

magnesium  and  1.493  in  the  case  of  stainless  steel.  Tokina  1.5  as  a 

representative  value  introduces  an  appro.* mat >11  enui valent  to  10  uncertainties 

in  the  dimension  loss  nower  per  unit  thickncs  ano 
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the  upper  solid  curve  of  Figure  2  is  found  to  result.  The  lower  solid 

curve  is  a  reproduction  of  the  results  given  as  Figure  1,  and  the  dashed 

$ 

curve  presents  a  smooth  interpolation  joining  the  two  solid  segments. 

The  resulting  composite  curve  is  suggested  as  a  means  of  rapidly 
estimating  the  time  reauired  for  complete  melt  through,  together  with  the 
definition  of  dimensionless  oower  per  unit  thickness  (Equation  3>3)  and 
the  one  dimensional  melt  through  time  (Equation  20). 

As  noted  earlier,  these  results  were  obtained  under  the  assumption 
that  no  heat  is  lost  to  the  surroundings.  If  heat  is  lost  from  each 
unit  area  by  connection  and  radiation,  then  the  temperature  must  satisfy 
the  differential  equation 


kV2T  spCoE  +  thihh)  + 

p  at  i  t 


(41) 


where  TQ  is  the  temperature  of  the  surroundings,  h  is  the  surface  heat 
taansfer  coefficient,  s  is  the  emissivity  and  o  is  the  Steohan-Boltzman 
constant.  An  approximate  solution  for  the  radially  sytranetric  steady  state 
temoerature  distribution  can  be  obtained  by  setting 

h*  *  h  +  4ea  T£ve  (42) 

which  yields 

T  a  T0  +  A  I  ( /2s  L  )  +  B  k0  (m  r-  )  (43) 

0  £  t 

L  ^  j 

where  8  =  —  ,  is  a  Biot  number  and  I.  and  K  are  modified  Bessel  functions 
k  u  r> 

of  the  first  and  second  kind.  For  small  arnumet.'s  the  following  approxima¬ 
tions  may  be  used: 
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l  ; 


i  f 


E  fc 


B  f* 


KQ(z)  «•  -  tn  z 


(44a) 

Kj(z)  -  1/z  (44b) 

Assuming  a  flux  FQ  is  again  uniformly  distributed  over  a  circle 
of  radios  R,  and  a  quasi  steady  temperature  to  exist  in  the  region  under  the 
spot,  then  a  heat  balance  on  material  under  the  spot  gives 

(45) 


F0irR2  =  2ttRi  (-k  —  ) 

0  3r  r=R 


Solving  for  the  constants  A  and  B,  we  find  A  =  0  as  the  temperatures  must 
remain  finite  at  large  r  and 


B 


FnR 


2k/ZH  Kj(/28  R ft) 

The  steady  state  temperature  distribution,  asymptotically  aonroached  at 
large  time  is 


(46) 


F0R  KJ/U  R ft)  „  _  €  . 

T  e  — ~ — °  _  ■■ — -  +  T0  for  r  >  R  ano  T  =  T(R)  for  r  <  R. 

^  I.  MV  1/  /  Ma  m \  ^  •“ 


2k/2{f  Kj(/2fj  R/i) 


(47) 


For  /2&  R/i  «  1,  the  asymptotic  values  may  be  used  to  obtain 


T(R)  ,  M_  ln  J  J__? 

i  U/2b( 


or 


(48) 


T(R)  =  T  ♦  It  Tm-Tc  +  id  1 


fi 

2s 


'P  J 


in  ? 


/  I 


\  R/2jf  I 


(49) 


The  critical  power  per  unit  thickness  required  to  bring  the  area  under 
a  uniform  beam  to  the  melting  temperature  is  therefore 


i  crit 


2s 


1 


1  + 


(50) 


Cn(T  -T) 
P'  m  o 


in  t/R/2? 
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An  identical  result  is  obtained  for  a  oaussian  beam  of  radius 

R  =  2a.  At  lower  values  of  power  per  unit  thickness,  melting 

will  not  be  achieved.  For  values  not  substantially  above  the  critical  value, 

the  times  predicted  by  the  method  described  in  the  previous  section 

(or  any  other  method  which  ignores  losses)  will  be  much  increased. 


v'a,  -  '.•Ait 


IV.  RETENTION  AND  VAPORIZATION  OF  MELT 


In  the  preceding  sections,  the  time  required  for  complete  melt- 


through  under  a  Gaussian  distribution  was  considered  and  compared  with 


the  times  which  would  be  required  assuming  one  dimensional  heat  conduc¬ 


tion.  The  assumption  of  complete  removal  of  the  melted  material  may  be 


justified  under  conditions  of  strong  air  flow,  but  the  alternative,  i.e.. 


molten  material  retained  as  a  liquid  until  the  vaporization  temperature 


is  reached,  should  also  be  considered. 


For  this  study,  the  finite  element  computer  program,  as  described 


earlier,  was  modified.  The  specific  heat  was  taken  to  be  constant  below 


T^,  the  temperature  below  which  the  material  is  solid,  and  above  Ty,  the 
temperature  at  which  the  material  becomes  liquid.  The  first  phase  change 
(melting)  was  accounted  for  by  assuming  an  effective  specific  heat  C 


over  the  m  lting  range,  i.e.,  between  the  lower  and  upper  melting  tempera¬ 


tures.  The  heat  of  fusion  is  assumed  to  be  added  uniformly  over  the 


melting  range  rather  than  at  a  fixed,  single  melting  temperature.  Let 


Cp*Cpo  for  T<TLorT>Ty 


C  »  C  +  L  /(T..-T  )  for  T,  <  T  <T„ 
p  po  m  U  L  L  U 


where  Lm  is  the  heat  of  fusion.  The  second  phase  change,  vaporization  at 


a  temperature  T^was  handled  as  melting  was  treated  in  the  original 
program,  i.e.,  the  finite  increment  of  mass  is  assumed  to  remain  at  Ty 
until  the  entire  heat  of  vaporization  is  absorbed.  At  this  time,  the 


entire  cell  is  assumed  to  be  instantaneously  removed. 


Table  III 

PHASE  TRANSFORMATION  PROPERTIES  FOR  SEVERAL  MATERIALS 


Material 

*L 

Ty 

TV 

Lv 

°K 

°K 

°K 

Ooules/g,, 

Aluminum 

916 

930 

2873 

'0492 

Magnesium 

878 

905 

1330 

5852 

Stainless 

Steel 

1673 

1728 

3273 

6360 

Titanium 

180’ 

1900 

Melt-through  tises  were  cos^uted  for  the  fear  z&tcris is  Tifcaniu*, 
Aluainua,  Stainless  Steel  and  Magnesian-  The  Gaussian  teas  was  assvaed 
to  have  a  diaeater  (4a)  of  1.04  cn  end  the  ssa»?e  thickness  was  assumed 
to  he  .03  ca.  For  Magr>esiuac  Stainless  Steel  and  Tttaniwn,  eh*  density, 
specific  heat,  conductivity  aed  heat  of  fusion  were  taken  as  given 
oreviously  (Table  i),  while  the  properties  of  aluninua  were  taken  to 
be  p  *  2.7  ge/ar*  K  *  2.01  J oule/oi  sec  *£  C.  -  .93G5  Ja#jle/ga*C 
l,.  *  395.4  Joules/er;.  Other  thervusl  properties  were  taken  as  given  in 
Table  III.  Coopyted  celt-through  tines  (t,,,  the  tine  at  which  the  rear 
surface  reaches  the  upper  welt  temperature  T #)  were  computed  and  are 
compared  with  tj,  the  tine  which  would  be  required  for  one-dtaensirna! 
heating  with  complete  removal  si  the  sane  intensity-  =h*  results  are 
tabulated  in  Table  I¥  and  shqwn  graphically  in  Figure  3,  ids ere  the  ratio 
of  celt- through  tiaes  is  given  as  a  function  of  rfieensferfess  p©*er  per 
unit  thickness.  The  solid  line  is  the  average  of  the  results  for  all 
materials  given  in  Figure  1,  and  represents  celling  with  instantaneous 
celt  renewal .  Fro*  Figure  3  it  can  be  see  ■*.  that  the  retained  melt 
causes  only  a  slight  increase  in  nelt-tfcrttfgh  tine  until  a  critical  valor 
of  dicensionless  powr  per  unit  thickness  >s  reached-  It  is  interesting 
to  note  that  this  value  is  comparable  tc  that  required  for  one-diner, sional 
heating.  At  larger  values  of  P.a,  the  tins  *-atfo  increases,  as  wre  and 
score  cf  the  incident  energy  is  "wasted"  in  evaporation  of  the  celt.  In 
this  case,  the  results  for  all  caterials  do  not  f«H  on  a  single  cure*. 
This  car,  be  attributed  to  the  fact  that  the  ratio  of  heat  of  vaporization 
to  heat  of  fusion  is  not  the  sans  for  z1 i  raterials  (Table  III).  The  re¬ 
sults  in  Figure  3  can  ba  seen  to  bs  ordered  as  the  ratio  of  beabiu 

?0 


Table  IV 


COMPUTED  MELT  THROUGH  TIMES  FOR  SEVERAL  MATERIALS 


Material  I  -Kw/cm^  tm-sec 

pa 

Titanium  .5  2.122 

•8  l!o69 

1  .807 

2  .3755 

5  .160 

10  .0965 

20  .0750 

30  .0640 

SO  .0505 

100  .0318 

Magnesium  1  > 2833 

2  ‘.0996 

S  .0355 

10  .0184 

20  .0112 

30  .0096 

SO  .0081 

100  .0059 

Stainless  Steel  1  .8688 

2  .3736 

5  .149 

10  .0820 

20  .056 

30  .0496 

50  .0423 

100  .0302 

200  .0192 

Aluminum  2  .228 

3  ’  ! 1166 

5  .0591 

10  .0273 

20  .0140 

50  .00622 

100  .0035 

200  .0024 

500  .0019 

800  .0017 


tj-sec 

1.151 

.7192 

.5754 

.2877 

.1151 

.05754 

.02877 

.01918 

.01151 

.005754 

.1373 

.06863 

.02745 

.01373 

.006863 

.004575 

.002745 

.00137 

.5629 

.2814 

.1126 

.05629 

.02814 

.01876 

.01126 

.005629 

.002814 

.1062 

.07081 

.04249 

.02124 

.01062 

.004249 

.002124 

.001062 

.0004249 

.000266 
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These  results  suggest  that  vaporization  effects  in  a  fully  retained 
melt  first  become  significant  at  a  value  of  power  per  unit  thickness 
which  is  comparable  to  that  required  to  achie  e  one-dimensional  axial 
heat  flow.  Consequently,  it  is  to  be  expected  that  tb)  heat  conduction 
aspects  of  laser  material  interaction  problems  nay  b?  treated  as  one 
dimensional  in  the  regime  where  vaporization  Is  significant.  Conversely 
it  <s  to  be  expected  that  the  vaporization  of  retained  melt  will  be 
negligible  in  cases  where  radial  conduction  is  significant. 


i 
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*  V.  MELTING  AND  VAPORIZATION  DUE  TO  A  SINGLE  PULSE 
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The  computer  program  as  modified  to  n»odel  two  dimensional 
heating  and  melting  with  the  melt  being  retained  throuqh  to  complete 
vaporization  was  further  modified  to  treat  single  or  repeated  pulses. 

For  pulses  long  compared  to  the  time  steps  used  in  the  finite  element 
method,  the  method  is  applicable  during  the  duration  of  a  oulse  which 
is  constant  or  which  increases  in  intensity  with  time.  In  order  to 
analyze  the  heat  transfer  and  melting  which  occurs  after  the  incident 
energy  is  set  to  zero  or  reduced  significantly  (as  at  the  end  of  a 
pulse),  account  must  be  taken  of  the  fact  that  the  temperature  of 
certain  cells  will  decrease  as  time  progresses.  A  cell  which,  at  the 
cessation  of  heat  addition,  was  in  the  process  of  vaporizing  must  either 
be  assumed  to  "give  up"  as  heat  flowing  into  other  cells  that  portion  of 
the  heat  of  vaporization  which  had  been  supplied  cr  else  must  be  assumed 
to  partially  vaporize.  The  latter  aoDroach  was  found  to  lead  to  smoother 
and  more  consistent  results  in  the  cases  considered,  particularly  when  the 
thickness  was  divided  into  only  ten  layers.  A  nulse  with  a  lenn,  slcwly 
diminishing  "tail"  would  reouire  special  treatment,  as  the  material  under 
the  beam  may  lose  energy  by  conduction  faster  than  supplied  by  the 
incident  radiation.  Repeated  pulses  present  no  difficulty  once  the  program 
is  modified  to  permit  "cooling  down"  between  pulses.  In  all  cases  considered 
here,  the  pulse  was  assumed  to  be  of  uniform  intensity  for  the  duration 
of  the  single  pulse. 

Limited  calculations  were  performed  ? 
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for  .08  cm  thick  sheets  of  aluminum  and  titanium,  assuming  gaussian 
beam  profiles  of  o  *  .179  and  .228  cm  and  material  properties  as  given 
in  the  previous  section.  Various  absorbed  intensities  were  assumed  in 
temporally  uniform  pulse  of  .005  seconds  duration.  Twenty  divisions 
through  the  thickness  were  used  in  all  cases. 

The  results  obtained  for  aluminum  showed  that  no  vaporization  is 
to  be  expected.  The  peak  absorbed  intensity  reouired  to  produce  melting 
to  various  depths  is  shown  in  Fig.  4.  Little  dependence  on  soot  size 
was  found,  as  is  evidenced  by  the  results  from  a  one  dimensional  calculation 
shown  for  comparison. 

In  the  case  of  Titanium  (Figure  5)  the  sheet  is  found  to  be 
partially  melted  and  partially  vaporized  by  the  beam.  The  portion  melted 
is  seen  to  increase  rapidly  with  intensity  until  the  threshold  for 
vaporization  occurs.  Above  this  intensity,  a  fairly  uniform  melt  thickness 
(about  .008  cm  in  this  case)  is  found  in  front  of  the  vaporization  depth. 

These  limited  results  suggest  that  the  spot  size  is  not  critical,  particularly 
in  regard  to  the  prediction  of  the  vaporization  depth. 

These  calculations  were  compared  with  experimental  data  obtained  from 
the  Air  Force  Weapons  Laboratory  in  order  to  determine  an  average  value 
for  the  absorbtivity.  If  a  gaussian  profile  is  assumed,  the  comparison 
suggests  an  absorbtivity  of  6.9%  for  the  experiment  on  aluminum  in  vacuum 
and  values  of  13.7%  and  9.2%  respectively  for  the  two  experiments  on 
Titanium  in  air.  The  beam  shape  in  these  experiments  would,  however,  have 
a  pronounced  effect  on  the  results  of  this  comparison.  Absorbtivities  of 
twice  there  values  would  have  been  necessary  if  th-e  beam  was  "flat,1'  while 
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the  observed  depths  of  melting  and  vaporization  would  be  possible  with 
a  lower  absorbtivity  if  the  vaporization  and  melting  occured  at  a  local 
"hot-spot"  in  the  beam.  These  limited  results  do  suggest ,  however,  a 
method  for  determining  absorbtivity  if  the  oeam  profile  is  well  charac¬ 
terized. 
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VI.  CONCLUSIONS 


Further  theoretical  studies  of  laser  induced  melting  and  vaporization 
using  the  finite  element  numerical  method  previously  described  have 
shown  the  method  to  be  adaptable  to  several  new  problems. 

Principal  new  results  are  as  follows: 

(1)  A  single  curve  has  been  developed  which  permits  the  estimation 
of  time  required  to  produce  melt-through  by  a  continuous  beam  of  gaussian 
profile  even  for  cases  where  radial  conduction  is  substantial.  The  melt- 
through  time  was  shown  to  be  a  unique  multiple  of  the  time  required  for 
one  dimensional  melting  at  the  same  intensity,  where  the  multiple  is  a 
function  only  of  a  certain  dimensionless  power  per  unit  thickness, 
independent  of  intensity,  spot  siz'.  and  material. 

(2)  A  simple  analytical  result  was  developed  which  shows  that  for 
values  of  dimensionless  power  per  unit  thickness  below  a  critical  value 
(which  does  depend  on  material,  spot  size  and  thickness)  melting  of  thin 
sheets  cannot  be  produced. 

(3)At  high  values  of  dimensionless  power  per  unit  thickness,  all 

problems  may  be  treated  as  0*0  dimensional  (axial)  heat  conduction.  The 

critical  value  is  about  P  =70. 

ta 

(4)  It  is  particularly  significant  that  the  geometrical  condition 
necessary  for  one  dimensionality  is  not  a  ratio  of  spot  dimension  to 
thickness,  but  rather  a  ratio  of  spot  area  (dimension  squared)  to  thick¬ 
ness. 

(5)  The  value  of  power  per  unit  thickness  necessary  to  produce  one 
dimensional  axial  heat  conduction  is  also  approximately  the  threshold 
where  vaporization  effects  first  become  significant.  At  higher  values 
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of  PAa  the  dimensionless  power  per  unit  thickness*  large  amounts  of 
input  energy  serve  only  to  produce  vaporization  of  the  melt. 

(6)  The  method  was  also  found  applicable  to  the  prediction  of  t.‘,e  deqree  of 
melting  and  vaporization  due  to  single  pulses  of  high  intensity  and 
provides  a  method  for  estimating  absorbtivities  if  the  melt  depth  due  to 
a  well  characterized  pulse  is  known. 
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